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The "bending stresses in the covers of "box oeans or 
de— flange Deans differ appreciably from the stresses pre- 
dicted "by the ordinary "bending theory on account of shear 
deformation of the flanges. The problem of predicting' 
these differences has "become known as the shear-lag problem. 

The first part of the paper deals with methods of 
shear -lag analysis suitable for practical use. The basic 
elements of these methods have been published in previous 
papers, but the treatment of these methods presented in 
this paper is consolidated and improved in several respects. 
The methods are sufficiently general to cover any arbitrary 
spanwise variation of cross section and loading as - oil as 
chordvise variations of stringer area, stringer spacing, 
and sheet thickness. Methods of analyzing the effects of 
cut-outs are also given'. 

The second part of the paper describes strain-gage 
tests made by the 1TACA to verify the theory. Three tests 
were made on axially loaded panels of variable cross sec- 
tion, six were made on beams of variable cross section, 
and three were made on beams of constant cross section for 
extreme or limiting cases. Three tests published by other 
investigators are also analyzed by the proposed method. 

In order to make the test of the theory as severe as 
possible, the NACA specimens were designed to show larger 
shear-la s effects than may be expected in typical present- 
day construction. The agreement was quite satisfactory 
even in extr.eme cases such as ver^ short, wide beams. Sat- 
isfactory agreement was also found in tests on the limit- 
ing case of a cover without stiff eners; this agreement 
shows that the theory is applicable to the case of heavy 
cover plates used without stiffening or to cases in which 
continuous stiffening in the form of corrugated sheet is 
used. 

The third part of the paper gives numerical examples 
illustrating the methods of anal-sis. An appendix gives 
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comparisons with |ther methods, particularly with the meth- 
od of Ebner and Koller. 



i::j^cruGTiGi T 



The "bending stresses in Vox be aft g do not always con- 
form very closely tc the predictions of the en r ;inocrin" 
theory of oondin^. The deviations iron the t'aeory are 
caused chiefly fcy the shear deformations in the cover of the 
oox that constitutes the flange of the bean. The problem 
of analyzing these deviations from the engineering theory 
of "oendin-; has he core knovrn as the shear-la? problem, a 
tern that is convenient although not very descriptive. 

The most important case of shear-la^ action occurs in 
the vr±&% structure. The cross section of the uin^ usually 
varies consideraoly alon~ the span; analytical solutions 
"based on the assumption of constant cross section are 
therefore of- little practical value, and methods of analysis 
have had to oc developed to cope with the conditions found 
in actual structures. The development of such methods has 
boon continued over a period of several ~ears (references 1 
to 3) and it is now possible to Rjivo a reasonably well- 
rounded presentation of practical methods of analysis. 

The paper is divided into three parts. The first part 
discusses the methods of analysis. The second part de- 
scribes tests made by the 1TA3A and shows comparisons be- 
tween experimental .and calculated results for the NAOA 
tests as well as for tests made elsewhere. numerical ex- 
amples to illustrate the methods of analysis are presented 
in the third part. 

The method of presentation chosen is intended to meet 
the needs of the practicing stress analyst. The paper con- 
tains the information actually needed in stress analysis; 
detailed derivations and discussions have been omitted, 
but they may be found in several of the cited references. 
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I . IS. E T I 0 D S 0 J A 17 A L Y S I S 
I>3? I III II 017 OP I HE PRCBLBI.' AND BASIC ASSUMPTIONS 

Seduced to itn simplest form the -orobler. r.ay be stated 
as follows: A sheet, stiffened or unstif f cr.ed, is fastened 
ton. foundation along one edge and loaded along the two 
edges perpenli cular to the foundation by distributed or 
concentrated forces as indicated in figure 1. The sheet 
say do a structure in itself (fig. 2(a)) or it v.. ay bo the 
cover of a box bear. (fig. 2(b)). The pro 1)1 or. is to find the 
stresses in the sheet. 

As shown in figure 1, stiff oners arc theoroti cally 
necessary along the loaded edges if concentrated forces ? 
are introduced because the stresses would otherwise becor.e 
infinite. These cd^e stiffencrs will bo referred to 
throughout this paper as "corner flaagos" or sir.aly "flanges." 
Otner stiffeners parallel to the loaded edges will be re- 
i erred to as " longitudinals" or "stringers"; these stiff on- 
ers nay or .nay n ot exist in aay & iren case and »a» 0 r nay 
not be attached to the foundation. 

It- will be assumed that the structure is always sym- 
metrical about a longitudinal plaao ( Tr = 0). This as- 
sumption materially si mpli f i o s the problor. without decreas- 
ing the practical usefulness of the theory very much be- 
cause r.ost practical structures are at least ao- roxinat &ly 
symmetrical. On account of the syr-.otry, it will bo effi- 
cient to consider one-half of the structure in all deriva- 
tion s and comput at ions. 

It will be assumed that infinitely many ribs of infi- 
nite exter.sior.al (chordwise) stiffness are distributed 
along the span. An equivalent assumption is frequently 
maae in theoretical solutions of stress problems. The t 
sumption is plausible in this case because it is fairlv 
oovious that the extensional stiffness of the ribs together 
witn the lateral ..bending stiffness of the flanges between 
the nos is sufficient tr take caro of such transverse 
stresses as might arise from longitudinal forces and stress- 
Cs, 1 T " e fir -al proof that the assumption of rigid ribs is 
admissible must, of course, bo furnished ** exoe rir.ent s 
like tnose described in the second port of this paper. 

The field of shear-lag analysis is very extensive; it 
was therefore considered advisable to confine the discus- 
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sion, in general , to "beans with flat covers. The no at gen- 
eral netlxod of analysis ^ivon in this paper can he rer? 
readily extended to heans with camhered covers and this 
extension is therefore ^iven. An approximate method for 
dealing with moderate anount s of canter is -^ivon in refer- 
ence 2. 

ANALYSIS OF SIK&LE~3TRIH&EH STRUCTURES 



Structures like those shown in figure 2, having hut a 
single stringer, are rarely encountered in iDractice. Nev- 
ertheless, the analysis of single-it ringer structures will 
he fully discussed for several reasons* The innediat e 
reason is that the fundamental relations as well as all 
the r.ethods of analysis can he easily demons t rat e d on this 
type of structure. A nore important reason is the fact 
that the most rapid method of analyzing multi stringer 
structures is based on the temporary reduction of the 
nulti stringer structure to a simple— strin^or structure, 

Sim Conventions 

The si ^n conventions adopted are as follows • Normal 
stresses and strains in the stringers and the flanges are 
positive when they are tensile. Shear stresses and strains 
in the cover sheet are positive when they are caused oy 
positive strains in the flange. Shear stresses in the web 
are positive when they are causing positive strains in the 
flange . 

The compression side of the "beam is analyzed independ- 
ently of the tension side* It is therefore permissihle and 
convenient to retain the sisn convention just ^iven for the 
analysis of the compression side, chan^in^ only the defini- 
tion of stringer stresses to positive when compressive. 

In general, the positive directions of the coordinate 
axes trill he taken as shown in figure 3. In some cases, 
particularly for analytical solutions, it is more conven- 
ient to use the opposite direction for the positive in- 
direction "because the resulting formulas are simpler.. 
(See, for instance, formulas for axially loaded panels, 
references 1 and 2.) 
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Fundamental Equations and Analytical Solutions 

For purposes of shear-la 1 ? analysis, all structures 
are idealised in a manner familiar, for instance, from the 
design of plate girders. Stringers arc assumed to be con- 
centrated at their controids; the idealized shoot is as- 
sumed to carry only shear, out the fact that the actual 
sheet carried lon^i tudinal stresses in addition to the 
shear is taken into account by adding the \70ll— kno-rm effec- 
tive v/idth of the sheet to the stringers. The participa- 
tion of the shear v;ob in the bending action is expressed 
"by adding l/o ht^ to A ? , 'rhich makes the section mod- 

iuus of the idealized section equal to that of the actual 
section. Figure 4 shows the idealized cross sections of a 
single-stringer panel and of a sin°;le-st ringer "beam; the 
standard "basic synods used in this paper are indicated in 
this figure. A complete list of symbols is -,'iven in ap- 
pendix A. 

Figure 5 shows an idealized sin^lc-st ringer bean of 
constant cross section subjected to a transverse load at 
the tip. Inspection of the free-body diagrams in figure 
5(b) shows that there arc two equations of static equilib- 
rium, 

m* = Stt ~ aSn . (la) 

c\F L = dS c (lb) 

where Srj is the shear force in the v:ob, in this case 
equal to P; and dS c = T t dx f where T denotes the 
shear stress in the cover shoot. 

Under the assumption of infinite transverse stiffness, 
the relative longitudinal displacement (up - u L ) of two 

corresponding points on the flange and on the longitudinal 
divided by the width b defines the shear strain Y and 
therefore the shear stress T (fi~. 5(c)). Because the 
displacements u are ^iven by the expression 

x=x 

r 5 * 

J jj dx, differentiation -"yives the basic elastic rela- 
x=L 
t i 0 n 

&t = - - G - (j - <j ) dx (lo) 
Eb * L 



whore Gr is the effective shear modulus, irhich takes into 
account the effects of "buckling -/hen necessary. Equations 
(la), (To), and (ic) can be conoinod to form a differen- 
tial equation, and this equation can do solved for simple 
cases. A number of solutions are s'iven in references 1 and 
2; similar solutions have "been : ;ivcn "by other authors. 
These analytical solutions are of some value in making 
comparative studies and in studying various aspects of the 
shear-la% proolcn. For practical stress analysis, hov/evor, 
numerical methods capable of dealing with arbitral^ vari- 
ations of cross section and loading aro required. Two such 
methods villi "be described: the solution by means of a re- 
currence formula and the solution "by successive shear- 
fault reduction 9 

Analysi s of Singl e-S t r Luge r 
Structures b7 r the Recurrence Formula 

ilinci^lQ and scope of ..method.- She principle of an- 
alyzing a beam of variable cross section is as follows: 
The "beam is divided into a convenient number of bays in 
such a ray that the cross section and the running shear in 
the web S^/h may bo assumed to bo constant within each 
bay. The shear deformation in the cover shoot of each bay 
is computed in terms of the unknown forces acting between 
bays. Application of the principle of consistent deforma- 
tions then : ;ivcs a sot of equations, similar in form to 
throe— moment equations, for the unknown forces. 

Theoretically, the method permits taking into account 
any variation of cross section and loading alon^ the span.- 
The limitations are t similar to those encountered in other 
problems of stress distribution in cases of variable cross 
section and loading. 

lftgtt3rrea.ee. formula f o zhfra? lag«»* As stated in the 
preceding section,, the beam is divided into a number of 
bays; tile cross section and the web shear Srj/h are as- 
sumed to be constant within each bay. The lengths of the 
bays need not* be equal nor need they be small, as is often 
required in similar methods. In the limit, a single bay 
may span the entire length of the beam. The system of num- 
bering the stations and the feasrs between stations is shown 
in figure 6. 

Each individual ba* r can nov; be treated as a free body 
subjected to certain forces (fi^. ?)• These forces can be 
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split into tv:o %T<mp® (£'!£« S): one ^roup consists of tho 
forces calculated by the ordinary bending theory, rhich 
assiiisss no shear def o mat ion ; tho other °;roup represents 
tho l:.:'ieronccs ootroo:: the actual forces an d the forces 
of tne first ^roup or ; in other v/ords, the changes in 
forces caused by the shear deformation of the cover sheet. 

The first ^roup of forces \vill bo designated P~forccs 
to indicate that they are calculated by the theory that 
assumes plane sections to remain plane. Individual forces 
and stresses ' belon^in^; to this ^roup vrill be denoted by a 
superscript P. The calculation of* these forces and 
stresses is familiar to every engineer and consequently 
need not oe discussed in detail. 

The second ^roup of forces will be designated X-forces, 
Because the P-forccs on any one bay are in static equilib- 
rium, the X-forces on any one bay must be a s elf -e qui li brat ed 
;-roup longitudinally; that is, at any %iven station the 
force X F acting on the flange must be equal and opposite 
to the force X L acting on the icn -;i tudinal . Thi s <>0EL^i^ .1^ 

sion v;as anticipated in figure 8 by writing X vri t hout ^ ' 
and L. 

The shear deformation of the cover sheet can nov; be 
calculated in terms of the known P-forcos and the unknown 
X-forces; the details of this calculation are -\'iven in ref- 
erence 3. ^ Equating the deformations at the adjoinin-; ends 
of successive bays yields the recurrence formula 

X n-i ^n " X P>n + Pn+i } + * n + 1 %H = ~ 7 n + (« ) 

where 

p = 2 

^ n Gt tanh KL 



Gt sinh KL 



(3a) 
C*T>) 



K I® - nhear-la^ par-meter appearing in all analytical so- 
lutions for sin^lc-st ringer structures (references 1 and 
2) and is defined bv 



X 



E"b 



(4) 



la equations (3a) to (3c) , each individual quantity should 
"be understood tc have a suhscript n, indicating the av- 
erage value for the hay in question. ilote should he taken 
that this statement applies to L, vrhich is to he taken 
as the length of the individual hay in question, not as 
the length of the entire "beam, 

Strictlv speaking, all coefficients Y appearing in 
this paper should have a superscript ?. These super- 
scripts have "been omitted "because they arc not needed in 
the actual use of the equations; they are needed only in 
the derivation of the equations (reference 3 ) • 

TTritten in more explicit form, the equations are 
Xo^i ^ x i(Pi * Ps ) + * + 7 2 

- x 2 (p 2 + p 3 ) + x 3 q 3 = + y 3 

. •.. . • ♦ i l • 

X n-^n ~ X n^n + P»+j) + X n-fi^i T ~ 7 n + " Vm ) < 5 > 



Z r-.^r - X r ( ^r + Pr+i> = ~ Y r + Y r+ i 

/ 

It v:ill he noted that the externally applied load appears 
only in the coefficients Y; for any ^iven heam, .then, 
the loft-hand side of the equations remains unchanged if 
changes occur in the loading. 

Boundary conditions . - Before the system of equations 
(5) can he solved, the "boundary conditions at the tip and 
at the root must he defined. At the tip, the following 
cases may arise: 

(1) Only a transverse force is applied (fi-$« 9(a)). 

In this case, Z 0 = 0. 

(2) A longitudinal force F may he introduced {fk%* 

9(h)). In this case . 

Xo = F±i (5) 

1 

0 
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When the longitudinal force P is the only force applied 
to the "beam, the idealized shear web is inactive, and the 
problem i g that of an axially loaded panel. 

At the root, the following cases may arise: 

1. The flan°;e and the longitudinal are connected 
to a ri^id foundation. 

i# The flange and the longitudinal are connected 
to a foundation that deforms under load. 

3* The flange is connected to the foundation; the 
longitudinal is not connected. 

^This system of classifying the possible cases is based 
on the convention of defining the foundation as the station 
where the vertical shear is taken out. 

Case 1 at the root arises in practice when a win-^ is 
continuous from tip to tip. The plane of symmetry i s 
equivalent to a ri^id foundation. Case 2 arises in prac- 
tice when a wln% is joined to car ry~throu^h members pass- 
i&q through the fuselage. Case 3 has boon used in prac- 
tical design to facilitrte the assembly of the to the 
fuselage by reducing the number of bolts" to a minimum. 

The foundation may be considered as bay r+1. In 
case 1 there is no shear deformation of bay r+1 , and 
p r+1 as well as Y r+1 equal zero. In case 2, Y r+1 
equals zero, because no shear is carried in bay r+1; the 
deformation of the bay depends only on the axial stiff- 
nesses of the flange and the longitudinal passing through 
the fuselage, and 

, C: ' . . ' : : : tn 

where L is the distance from the w&n% root -to the Diane 
of symmetry of the airplane. 

In case 3 the last equation of the system cannot be 
used, and X r is found by inspection to be 

M At 

*r = r- T (8) 
h T7 A T 
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Cal culat i o n_ o^strez s e s_ f ro.m „X r f o r c e s . - After t he aye-* 
t em of equations (5 ) has been solved, the lon^ituclin-al 
stresses -re. found &y superposing on the stresses ealeulat- 
ed b* r the ordinary bea&in^ formula the stresses calculated 
from the X—f orces 



and 



a F = (J* * X/Aji . (9^3 
d L - ft * - X/A L (9b) 



where is the stress calculated by t>.e ordinary hendin^ 

formula. In the case under discu^ion, v:horo the beam has 
no cam 'bar § 

. ^ . JL- do) 

The running shear in the cover sheet of "bay n close to 
the inboard end of the hay, that is, close to station n, 
is siven by the formula 



(xt) : rJ±L \ + 1 , — 5a x — *a — <iia) 



Hear the outboard end of bay a, that is, near station 
n - 1-,' the running shear in the cover sheet is 



<">». - (s%) + *»- t^srfe - x » s^fe (11 " 0) 

x n 

For some applications it is desired to compute the 
average running shear in a hay. If the hay is not too 
Ions, tMs .&vera?e shear ma? he obtained by avera^in^ the 
shears at the tv:o ends of the oay computed by the formulas 
just ^iven. The result is 



2„ K, 



< T V>„ - (s-^) - I « «. - W + «•> <»«> 

x n 



An alternative v:ay to compute the average shear is to use 
the basic static relation (lb) 
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(TtL) n = f» - Ft (lid) 

n— 1 

Formulc (lid) Sires the true average, while formula (lie) 
i s approximate • 

Influence of foppr Ip. A^-pth apA width**- When a "beam 
Is tapered in depth, it is necessary to remember that part 
of the vertical shear is carried Tgjr the inclined flanges 
and longitudinals, so that 

where i is the inclination of the tension flange \-:ith re- 
spect to the compression flange. 

When a beam is tapered in width, neither the ordinary 
□ending theory nor the shear-la? theory is strictly appli- 
cable. The error caused by applying the ordinary bending 
theory, however, is small for normal angles of taper; to a 
similar decree of approxinat i on , the following approximate 
method of shear-la? calculation nay be used. 

Assume that the taper is removed "by nakin? the widths 
b at all stations equal to the width $ at the root. 
At the sane time, increase the sheet thicknesses in the 
ratio b r /b. The result rill he an untapered beam that 
has the same shear stiffness Gt/b at any station as the 
actual beam* This method of procedure assumes that trans- 
verse components of longitudinal forces can be neglected; 
this assumption is in keepin? with the assumption of ri^id 
ribs. 

, I* should be noted that the parameter K (equation 

{4;})JLm any bay of the fictitious untapered bean is equal 
to the corresponding parameter X of the actual tapered 
beam, but the coefficients p 5 q, and Y of the ficti- 
tious beam differ from those of the actual beam by the 
ratio b/b r . It is stated in reference 3 that the effect 
of taper in plan form jai?Jit be more pronounced than is in- 
dicated by the method just ?iven. Re-exaninat ion of the 
test data in the li?ht of the additional test experience 
gained since reference 3 was written tends to show that 
the method -^iven here is sufficiently accurate for the ta- 
per ratios likely to be encountered on win's. 
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Analysis of Single-Stringer Structures ^7 
Successive Shear-Fault Reduct ion 

Principle ^ nd scope of 'method,-* The principle of t he 
method of successive shear-fault reduction is as follows . 

An estimate is made of the stresses in the flair 3 ; 

the stresses o" L in the longitudinal are calculated by 
statics. By the application of the oasic equation (1c) 
and a. process of numerical integration, the spanwise dis- 
tribution of shear force in the sheet can then he calcu- 
lated- On the other hand, application of the oasic equa- 
tion (lb) also sives a spanwise curve of shear force in 
the sheet. The tvro curves will not a^ree except by acci- 
dent "because the estimated values of dp and OV will 
not fulfill the elastic relations and the boundary condi- 
tions except by accident. The difference between the tvro 
curves trill be referred to as the curve of " shear fault s • 11 

The existence of shear faults in the calculation 
proves that the assumed stresses Op do not constitute 

the true solution of the stress problem for the specified 
exte-rnal loads. The assumed stresses Oji constitute, 

however, t&e true solution for a closely related problem, 
namely, the structure subjected to the specified external 
loads and, in addition, subjected to a system of external 
loads equal to the shear faults. Obviously, then, the de- 
sired solution can be obtained from the assumed solution 
by deducting the effects of the shear faults. This deduc- 
tion is effected by superposing the effects of corrective 
external shear forces that are assumed to be applied in 
opposite direction to the shear faults. 

If the magnitudes of the corrections were made equal 
to the faults, the basic static equation (lb) would be 
fulfilled at each station but the basic elastic relation 
(ic) ^.Tould be upset. As a compromise between these con- 
flicting requirements, the :orrccuon is made equal to 
one-half of the fault. 

3ecause transverse forces are absorbed by the rib 
system and are not considered, the introduction of an ex- 
ternal shear is equivalent to the introduction of a pair 
of equal a:;d opposite forces. By St. Tenant 1 s principle, 
the influence of such a combination of forces is felt over 
only a limited distance. In order to simplify the computa 
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felon, it -ill "bo assumed that tho influence of each cor- 
rective force decreases to zero at tho next station. Er- 
rors introduced by this simplification -ill he small and 
will eventually be eliminated by repeating the process of 
correction. 

Application of tho corrective forces to the Initial- 
ly assumed values of (T-c and ov yields a new set of 
values for c% and cr^ , and the entire process is re- 
peated. It will he found that tho corrective forces are 
oecomin 0 ; smaller with each repetition of the process, so 
that the solution will he obtained by a sufficient number 
of repetitions. In theory, the computation is finished 
when the corrections to Op and cr^ are reduced to one 
unit of the last significant figure of o* ? or o^. In 

practice, the computation will often he finished sooner at 
the discretion of the analyst. 

For sin-^le-st ringer structures, the method of succes- 
sive shear-fruit reduction is unlikely to he favored over 
the recurrence formula because the time required for a 
solution depends very much on the ability of the analyst 
to make a -;ocd initial estimate of cr-p and <r L . The time 

required for a solution by means of the recurrence formula, 
on the other hand, is almost independent of tho skill and 
the experience of the analyst because the only item left 
to his choice is the number of bays. The method of shear- 
fault reduction for sin^le-st ringer structiires, however, 
is the direct basis of the most general method for analyz- 
ing r.ultistrin?er structures, .and this . fact justifies the 
description of the method. 

Method of sue ce s alve shear*»£ault redtic t i^ n . - In order 
to apply the method of shear-fault reduction, the beam is 
divided into a convenient number of bays. Because the 
computation involves n\imerical integration and differenti- 
ation, the lengths hx of these bays must be chosen fair- 
ly small so that no appreciable error is made by assuming 
the stresses to vary linearlv in each bay. Five fcajrs may 
be considered as the minimum. In order to reduce the time 
required for computation and the possibility of errors, 
the bays should be made of equal lengths whenever feasible. 

Tho computation is started hy tabulating for each 
station the c >iven magnitudes of Ay, At, , t, S , and H/h 
(or P) if they vary alon^ the span. If the beam tapers in 
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width, a fictitious beam of constant width Is used, as pre- 
viously di scuss 3 d. 

The magnitudes just onunoratod should "bo separately 
tabulated oocauso they will remain constant'; whereas, the 
main part of the calculation is repeated a number of times « 
The details of the procedure are learned most easily by 
following column for column the numerical example ^iven in 
part III, table X* 

Column 1 in table X $1 70S assumed values for cr^ . In 

assuming these stress values, the analyst must be guided 
by previous experience. It is possible to use any arbi- 
trary values whatever out, if the assumed values differ too 
much from the true onos, a lar^e number of cycles of the 
computation will "be required. The simplest procedure for 
general use is to multiplv the stresses obtained from the 
ordinary bending theory oy a factor slightly larger than 
unity, tfith some experience, this factor can "be estimated 
reasonably well from a knowledge of the average of the 
shear-la^ parameter ICL and the loading condition. 

Column 2 ^ives the forces Fy = (FffAj . 

Column 3 ^ives the forces F* » M/h - Fjj in the case 
of a beam or Ft - ^ — 31j in the case of an axially load- 
ed panel. 

Column 4 ?ives the stresses a% == ^l/At,* 

Column 5 ^ives the differences between columns 1 and 

4 £e$ - a L ). 

Column 5 ^ives the increments of shear stress obtained 
from the basic relation (ic) , 

At = - ||S (Oj - a L ) (SS-i) 



It nrill be noted that the values of AT in column 6 are 
positive. This s± >B arises from the fact that the inte-. 
^ration of the shear stress increments proceeds from the 
root to the tip so that the increments Ax are negative. 

Column 7 %ives the shear stresses T in each bay. 
These stresses are obtained b?r adding up the increments 
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^'ivon in column 6, starti&g at the root where T = 0. It 
should "bo noted that the values of AT represent the In- 
crements of shear stress for intervals of length Ax alon^ 
the s p an ; the fi i s t an ce between the root &n d the ml d die of 
bay r is, however, cmly half an interval Ax, so that 
the value of T in the root bay is T = § A T. From here 
on, the full value of At is. added each tine, unless the 
value of t at the tip is to be calculated v:hon a one- 
half step would be used a°;ain. (The value of T at the 
tip is needed for the calculation of the margin of safety, 
but it is not needed for the calculations indicated in 
table X. Consequently, this value is calculated only after 
the last cycle has been completed.) 

Column 3 -^ives the increments of shear force 

AS C3 == Tt Ax (3S-2) 

Column 9 ^ives the increments AP~ , obtained by sub- 
tractin^ the value of Ft at the outboard end of the bay 
from the value of BV at the inboard end of the bay. 

According to the basic relation (ic ) , AI? L should 
equal ASq-™ in each bay. The difference^ in each bay con- 
stitute the shear faults 

S? = ASn T - AF. (SS-3) 

and the shear faults 31 are ^iven in column 10, 

Consider no\7 figure 10(a), \rhi ch shov/s one bay with a 
positive shear fault SI and the corresponding shear- 
fault correction SFC; S?C is in the form of external 
forces distributed uniformly alonn; the bay. 

The length of a bay is small compared v/itb the length 
of the structure; it may therefore be assumed that the 
properties of the structure just outboard and just inboard 
of the bay considered are the same. Under this assumption, 
one-half of the shear-fault correction SFC will be ab- 
sorbed by the structure outboard of the bay; the other 
half -..-ill bo absorbed T>7 the structure inboard of the hay. 
As previously stated, the total shear-fault corrective 
force '-.-ill be taken as one-half of the shear fault 



Total SFC » ~ t s ? 



(SS-4) 
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Therefore r-he corrective force at the outboard end of the 
"bay will "be 

SFC Q a <* | SFC a | SF (SS~5a) 

and the corrective force at the inooard end of the "bay 
v ill ce 

S?C i = |. SFC a * I SI 1 (SS-5h) 

The corrective stresses A dp and Acr^ are found by divid- 
ing the corrective forces SFC Q and SFC.- by the areas 

or 4t and are shown in figure 10(o). The si^ns of the 

corrective stress ACp ^re the ci c ;ns ^iven in formulas 

(SS~5a) and (SS-5"b) f while the Si?ns of the corrective 
stresses Ao^ ere opposite to those ^ivcii in /formulas 

(3S-5a) and (SS-5b) (fi». 10(h)). 

At the tip station there is no outboard structure to 
develop an' r resistance to the shear-fault correction force. 
Cons e cuentl y , for the tip "bay 

SFC Q S3 0 ( SS~5c) 

S7G ± a SFC = | 3? (SS~5d) 

In the numerical example, tahle X, it will "be seen 
that column 11 lists the values of SFC Q and column 12 

lists the values of SFC^ . At each station there is one 

value of SFG 0 and one value of S?C i . The sum of the 

two valties is the final value of the shear-fault corrective 
force and is tabulated in column 13. 

Column 14 fives A(7p - SFC/A-?, and column 15 c ;ives 
Aov = ~ SFC/A L . 

The addition of the corrections A Op to the initially 
assumed values of oy and of the corrections Ad^ to the 
initial values of OV ^ives a new set of values for Op 
and cr~. The cnti.ro process is then repeated r.s indicated 
in tabic XI jut the column °;ivin^ Fp is no longer needed* 

The entire calculation as shown in table XI is repealed 



17 



a^ain and a^ain until succossivo sots of values of Oj 
and c- -re judged to a-^ree ;/ith sufficient accuracy. 
The limit of possible accuracy Is reached v/hen the values 
of Ao-rp or A(J L hocome equal "to unity in the last sig- 
nificant figure of &f or ov. 

In or-der to avoid carr^in^ alor.-; errors, should 
"be obtained from the static equation F L m M/h - fj, every 
second or third cycle instead of from Aav . 

The sum of the shear faults mav "be used as an indica- 
tion that correct $i<*n conventions have been used; the 
sun of t : :.e faults in any -^iven cycle must "be smaller than 
the sun of the faults in the preceding cycle. This crite- 
rion is not sufficiently sensitive to prove the absence of 
any numerical error, out it is sometimes a v/el.-ome help 
\rhen starting calculations . 

A complication arises v/hen the longitudinal is not 
connected at the root. In this case, the stress <j l is 

equal to zero at ;;ho root out the shear stress t is not 
equal to zero. It is therefore impossible to proceed di- 
rectly v/ith the summation of the increments AT. In order 
to overcome this difficulty, a trial value T 0 for T at 
x = 0 in assumed, and the summation proceeds from this 
trial value. Prom statics, it is evident that 

x=L 

The trial value T q must therefore he negative, so that 
the summation of the increments AS fix aloBg the entire 
span ma? 1)0 equal to zero. On the first trial, this condi- 
tion trill not oo met except by accident, and the trial 
value for T q must be adjusted until the %irm condition 

is net. Speaking graphically, the process consists in 
finding the area be two op a curve (the Tt- curve) and an 
arbitrary horizontal line, and then shifting the horizon- 
tal line unr.il the area becomes zero. After the first ejr~ 
clc has boon compl-tcd, the value t q ootained can he 
used as a trial value for the second cycle, and it will 
oo so close that the necessary adjustment will he small. 

When the longitudinal is discontinuous at some x>oint 



18 



other than the root, the summit ion of the increments AT 



nay "be performed 
the root and the 
the outboard end 
m an n e r anal o % o u s 
di s continuous at 



in the usual manner for the region "between 
inboard end of the "break. The region from 
of the "break to the tip is treated in a 
to that just discussed for a longitudinal 
the root . 



In a cambered "bean, the "basic equation (lo) must be 
modified to r e a d 



dT = ~ 



Eb' L 



(cr, - o%) - (0, 



) 



dx 



a-; shov;n in reference 2. In this equation, Op is the 
stress in the flange calculated by the -usual Mc/l for- 
mula, and is the stress in the longitudinal calculat- 

ed by the Kc/l formula. In the case of a flat cover, 
oy equals G" L ? and they cancel , reducing equation (lc 1 ) 
to equation (lc). V/hen a "beam is analyzed "by the shear- 
fault-reduction method, formula (33-1 ) must he modified to 
conforin with formula (lc 1 )* An additional column will 
therefore be required after column 5 in table X. 
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ANALYSIS 0? MULT I S T EI KG- EH STRUCTURES 



Two methods will he -iven for the analysis of multi- 
stringer structures. The first method consists in reduc- 
ing the prohlem to that of a fictitious sin^le-strin^er 
structure that can he analyzed, "by the recurrence formula,. 
The final step of transferrin 0 ; back to the actual multi- 
stringer structure can he made onl rr under the assumption 
that - the chordwisc distribution of material - stringers 
and sheet - is uniform and that the moduli S and & are 
constant alon°; the chord. Small variations from uniform- 
ity can he disregarded cut , when lar^e variations exist, 
it is desirable to have a more general method available « 
For such cases, a method of successive shear-fault reduc- 
tion is described that is an extension of the method of 
successive shcar-f fruit reduction described for single- 
stringer structures. This method permits taking into ac- 
count arbitrary chordwise variations of stringer size, 
.stringer spacing, sheet thickness, and elastic moduli. 
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Substitute Simple— St ringer Method 

.?Zi^ciplG_q f .method.- The transverse "bending loads 
noting on a "box bo an are taken up first by the shear webs. 
The shear stresses in the w„eb are partly converted into 
normal stresses at the flan-^o; the rest of the stresses 
"becomes shear in the cover sheet., which is gradually con- 
verted into normal stresses in the longitudinals as the 
longitudinal plane of syinnet ry Is approached. It may he 
said, therefore, that the most important physical action- 
centers around the flange "because the conversion of shear 
stress into bending stress "begins here. 

'This consideration leads to a very convenient method 
of analyzing a mult i stringer structure by substituting 
temporarily a fictitious sin : ;le-strin^er structure. This 
fictitious structure retains without change those parts of 
the actual structure in which the primary and the most im- 
portant action takes place, namely, the shear web, the 
corner flange, and the sheet adjacent to it. The longitu- 
dinals, however, are combined into a single fictitious 
stringer, the "substitute single stringer," located at the 
centroid of the internal forces in the stringers. The anal 
ysis of the resulting sin ^le-strin~ er structure can be per- 
formed by the methods described previously and "ives the 
actual stress in. the flange (equation (9a)) as well as the 
actual shear stress in the cover sheet next to the flange 
(equations (ll)). For the stress in the longitudinals, 
only an average value is obtained by the analysis of the 
fictitious. sinTle-st ringer structure. The stresses in the 
individual longitudinals of the actual structure are calcu- 
lated at any ^iven station alon^ the span by assuming that 
the average stress just calculated is distributed chord- 
wise according to the hyperbolic-cosine law found in such 
analytical solutions as have been published. 

•The validity of the substitution method outlined can 
be made plausible in a general way by reference to St. 
Tenant ! s principle. A' :.:Uch ."ore convincing proof, how- 
ever, will be ^iven by the comparisons between experimen- 
tal and calculated results in the second, part of this paper 

^eterrii nation^ 
£]iie___(f irst ..approximation^ .- A typical cross section of a 
multi stringer structure is shown in figure 11(a). This 
cross section is idealized as indicated in figure 11(b). 
It should be noted that the effective width of skin adja- 
cent to the flange is considered as a longitudinal distinct 
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from the flange. The adoption of this rule makes it feasi- 
ble to cover all possible cases with a single rule because 
in a limiting case such as sho/.vn in figure 1, for instance, , 
obviousl/r the entire sheet should "be considered as consti- 
tuting the longitudinals. Incidentally, this rule tends 
to reduce the error due to the finite number of stringers 
that will be discussed. 

The width d of the idealized sheet between longitu- 
dinals depends on the spacing b x between rivet rows and 
on the type of the stiff oners. Open-section stiff eners 
(fi4f« 3,2 (a)) do not contribute to the shear stiffness of 
the cover; therefore, d = b x . Closed-section stiffeners 
(fi%. 12(b)) contribute to the shear stiffness of the cover. 
If this contribution is taken into account, the idealized 
Width for shear deformation is d = b x + b 2G , in v:hich. 

b O0 - ' (13) 

tp 

where t^ is the thickness of the stiff ener and p the 

perimeter, or developed v:idth, of the stiff ener "between 
rivet rows. 

The idealized mult i st ringer structure (fi<% 11(b)) is 
now converted into a siri^le-st ringer structure by combin- 
ing all idealized longitudinals into a single longitudinal 
located at the force centroid of the longitudinals. Be- 
cause the actual stresses are not known at this sta-;e, the 
stresses computed by the ordinary bending theory are used 
to obtain a first approximation. For the flat covers un- 
der consideration, the -'orcc centrcid will then be the 
centroid of the cross-sectional areas of the stringers, 
the Mc/l stress bein^ the same in all stringers. The 
distance of this centroid from the flange is the width bg 
of the substitute structure ('fig, 11(c)). The substitute 
structure can be analyze_d by the recurrence formula or by 
any other method if desired. If a second approximation 
is to be made, the calculations made for the first appro xi- / 
nation can be confined to finding the stresses op in the 
flange and cr^ in the longitudinal of the single- st ringer 
structure . 



the substitute single-stringer structure furnishes the 
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flange stress oj* and the chor&vrise average of the stress- 
es in the longitudinals for all stations along the span # 
The actual chordy/ise distribution of the stresses nay he 
obtained in the following manner, as explained in reference 



For the Uniting case of infinitely nanv stringers, 
some analytical solutions have been obtained in the form of 
solutions for the continuous cover sheet. These solutions 
shov: that the chorduise distribution of the stringer stress 
es at any given station follows a hyperbolic-cosine lav. 
The stress at a distance y from the center line nay there 
fore be vrritten as 



(J = 



°CL 00 sl - Y ^ 



(14) 



an auxiliary oaraneter and (jn r 



the value 



wnere i 

of o* at jr = 0. In this equation, both the stress c CL 
in the longitudinal at the center line and the auxiliary 
parameter Y are unknown 



Tvro conditions are availably 



to determine these unl^nov/ns: (l) The average of the 
stresses a between y = 0 and y = b nust be equal to 
the stress cr T of the substitute single stringer, and (2) 



at the flange y 



the stress 



nust ecual the stress 



Op. The result is a transcendental equation for Yb , 



Yb 



In order to facilitate the dot e rninat i o n of Yb , figure 
13 has been prepared. With the help of this figure, Yb 
can be determined by inspection after computing the ratio 
°Ij/ °]T * ^~' le s "fc res s a,t the center line is then computed by 
the formula 



(Jp/cosh Yb 



(15) 



In order to compute the stress in any stringer at a given 
distance 'y from the center line, it is only necessary 
to compute Yy = (Yb) x (y/b) and to apply formula (14). 

Formulas (14) to (lo) apply only whon .0 < a L / ' J? < 1. 

In regions critical for design rorlc, this condition is prob- 
abl: 5 * always fulfilled. For cprtain purposes such as check- 
ing the thcor:^ against experinental results, however, it 
ma'- be desired to calculate the chordvrise stress distribu- 



22 



tion net stations where the ratio o^/op falls outside of 
this range. It was proposed In reference 1 to replace for- 
mulas (14) to (16) for such cases "by 

cr = cr CL (2-.cosh Yy) (14a) 
sinh Yd 



Yd 



2 - cosh Yo °"y 



(15a) 



°CL = Sjt/ (2 - cosh Tb) (16a) 

Formula (15a) was used instead of formula (15) to extend 
the range of the Yh-curve in figure 13. It will he noted 
in figure 13 that the Yh-curve for very small negative 
values of C-^/ (J v does not become infinite as would he ex- 
pected by analogy with small positive values. This peculi- 
arity is caused oy the approximate nature of equation 
(14a) and is of no practical importance, 

' Qpr r egtion of c&ordwise stress di3tribuMon_f^^inile 
n^;D^^._o f_ stringer s . - The method of computing stringer 
stresses by using formula (14 ) is "based on the assumption 
that the stringers are. infinitely closely spaced. If the 
spacing of the stiff oners is finite, the total internal 
force -.-ill he found h. rr a summation instead of an integra- 
tion, and the internal force will differ somewhat from the 
external force. The iiagnitude of the error depends on the 
numher oi stringers and on the curvature of the chordwise 
stress plot, which is characterized "by the ratio (J^/cr-p or 

by the parameter Yh. 

The sign of the error depends on the location of the 
first stringer near the flange. Under the rules given for 
idealising the mult i st ringer cross section, the first full- 
size stringer is located at y = "b(l - l/n) , where n is 
the numher of st ringer s ( arrangement A, fig. 14). For 
this case the summation of the stringer forces will yield 
a smaller force than is necessary to halance the external 
load. If a full-size stringer were located at the edge 
y = h (arrangement 3. fig. 14), the summation of the 
stringer forces worild yield too large a value. As long as 
Yh is less than about 1#S , the errors for these two cases 
are numerically equal and are shown in figure 14. 
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The rule that the effective width of skin adjacent to 
the flange should oe considered as a stringer (fig. 11(d)) 
helps to reduce the error by Winging the actual case be- 
tween the txro oxtreme arrangements A and 3 of figure 14. 
In practice, the ratio of the actual force to the summation 
of the calculated stringer forces may be applied as a cor- 
rection factor to the calculated stringer stresses as illus- 
trated "by the numerical example in part III. This method of 
correction was used in the an-lysis of all NAG A tests de- 
scribed in part II with very satisfactory results, even in 
some quite extreme cases; it was also used with very satis- 
factory results in making comparisons with the Ebner-Koller 
method. (See appendix 3.) If the results obtained "by this 
method should he considered as too inaccurate, the method of 
successive shear-fault reduction may he resorted to for im- 
proving the accuracy of the results. 

Sue c e r .^i ve_aE2.?o xi mat ion b for _ su b s t i t ut e width . - By def- 
inition, the substitute widtb is the distance from the flange 
to the force controid of the stringers. For infinitely manv 
stringers, the controid can be found by integration (refer- 
ence 2), and its location is shown graphically in figure 15. 
The substitute width is ^iven by the expression 

In- any given case, the factor 1 - (y^/i ) is taken from 

figure 15. and b is the effective vidth for shear deforma- 
tion as defined by figure 11(b). 

In order to use figure 15 it is necessary to know the 
parameter Yb ; for this reason it is necessary to make suc- 
cessive approximations. In the first approximation it is 
assumed that there is no shear lag; in this case Yb = 0 and 
1 *■ '3 r L /o) = 0.5. The first approximation to the substitute 
width is therefore b s s f, and with this v/idth the first 

analysis is carried out as previously discussed. The stress- 
es C7p and (Jr are calculated for the substitute single- ( 
stringer structure, and for each station the ratio C" L /a F 
is calculated and used to determine the value of Yb from 
figure 14. The spanv:ise average of Yb is then calculated 
and the corresponding value of 1 - (y L /b) is found from 
figure 15. Tfci s nev; value of 1 - (jv/tl) is inserted in 
i ormula ;17) to obtain the second approximation to b<* , and 
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the analysis of the substitute sing;Ie>-st ringer structure 
is repeated with the changes necessitated bgr obliging the 
substitute width. 

If the stresses 0*^ and 0% obtained In the second 
approximation differ very ifiuch f row the stresses obtained 
in the ^irst approximation , a third approximation may "be 
made* On account of the rapid convergence of the process, 
the difference cet^ocn the first and the second approxima- 
tions need not "be very small to insure that the second 
approximation may be taken as final. It is suggested that 
the stress analyst work some examples "by means of the ana- 
lytical formulas given in reference 2. As a rough 6 ui do , 
it nay he stated that , if the accuracy of the 10-inch 
slide rule is used as a criterion, the second approxima- 
tion may oe considered as the final one when the shear-la^ 
parameter XL for the entire "beam is b rcater than 4 in 
the first approximation. 7Jhen Ki is about 7 or greater 
than, 7 in the first approximation, the f i rst approximation 
is sufficiently -accurate. These relations are also influ- 
enced to some extent "by the ratio A^/A^. 

The outlined procedure should "be r* lightly modified for 
axially loaded panels. In such panels, the value of Yo 
"becomes infinite at the station where the axial load is 
introduced. In order to avoid this difficulty, the span- 
wise average of the ratios cr- / °7 should "be found and Yb 
for the average ratio o>/ <Jj should he determined. This 

method may he applied to beams in many cases and the final 
results ohtained "by the two methods v:ill he the same, at 
least for practical purposes. It is preferable, however, 
to use the two distinct methods to avoid uncertainties in 
pro ce dure . 

The method given for flawing successive approximations 
to bg applies directly only when there are infinitely 
many stringers. When there are only a few stringers, the 
first approximation bg is not equal to b/2 out is de- 
termined oy the cent roid of the areas of the stringers as 
discussed in connection with figures 11(h) and 11(c). In 
such cases, it ma^ he assumed that the- ratio of a higher- 
order approximation of og to the first approximation 

hg^ is the same as though there were nan? stringers; any 

higher-order approximation to the substitute width is then 
given hy the expression 
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where the factor l-(y L /b) is determined as before from 
figure 15. 

Method of Successive Shear-Fault Reduction 

griftci^e &s& of pet&Q&»« The anal-sis 6f twit** 

stringer structures bv successive shear-fault reduction 
employs the sane "basic procedure that is used for the analy- 
sis of singlc-st ringer structures. Some modifications and 
additional concepts are, of course, required to adapt the 
method to the riuch more complicated problem of analyzing 
mult i stringer structures . 

The process of successive shear-fault reduction in a 
single-stringer structure consists in a repetition of ad- 
justments on a spanwise sequence of elements. It is obvi- 
ously not feasible to carry on such a process of adjust- 
ments at the same time on chordv7ise sequences of elements. 
In order to overcome this difficulty, a concept will bo 
introduced that has become quite familiar through the 
Cross method of moment di s t ri hut i on , namely, the concept 
of locking parts of the structure in place to isolate the 
part "being adjusted from the rest of the structure. The 
particular method of locking employed herein consists in 
locking certain stringers at a given state of longitudinal 
strain or, to use a descriptive expression, in imagining 
them to "be frozen solid. The stringers locked at any given 
time are the stringers to either side of the one "being ad- 
justed. The stringers are adjusted in sequence, starting 
from the flange and proceeding to the center-line stringer. 
The process is repeated until the agreement "between sue— 
ccssive cycles of the computation is considered sat i sf acto ry . 

The method is obviously more laborious than the sub- 
stitute single-stringer method. It is very general, how- 
ever, being capable of taking into account chordwise vari- 
ations of stringer spacing, stringer area,- sheet thickness, 
and shear modulus; it can also deal more successfully with 
structures having a very small number of stringers (two or 
three ) . 

In practice, it will probably be found advantageous, 
in general, to use the substitute single-stringer method 
to obtain a first approximation, average values being used 
wherever necessary. The method of shear-fault reduction 
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can then "bo used to improve the ?xcunc;- r of the results, 

5 lie method of shear-fault reduction has one advantage 
that na" be helpful at felines. After the constants have 
been computed and the first cycle has been completed, the 
work involved in sAicceeding cycles is so simple that it 
can he handled by computers with little engineering train- 
ing. 

Procedure, f or comp utation*-! The computation is started 
by assuming initial values for the stresses in all string- 
ers A to 3F (fig. 16), taking care that at each station 
the summation of the internal forces equals the external 
force U /h or P . 

The flange A is adjusted first. In order to effect 
this adjustment, the stringer 3 is locked at the state 
of stress initially assumed. The computation then proceeds 
in practically the sane manner as described for single- 
stringer structures; the only difference is that the val- 
ues of (in this case o'^ ) are not changed hut remain 
the same for all cycles. After a number of cycles - say 
five cycles - the adjustment of stringer A is stopped, 
and stringer A is locked at the state of stress just corn- 
put c d • 

Before the adjustment of stringer A was started, 
static equilibrium existed between the internal stringer 
stresses and the external load at each cross section. 
After the adjustment, equilibrium no longer exists; oof ore 
the adjustment of stringer 3 is started, it will he nec- 
essary to restore this equilibrium • To this end, the 
stresses in stringer 3 arc increased or decreased so that 
the summation of the internal forces at each station again 
equals the external force, 

With these corrected stresses acting in stringer 3, 
the adjustment of stringer 3 is started. Stringer . A ^is 
locked at the stresses obtained from the first process of 
adjustment; stringer 0 is locked at the stresses initial- 
ly assumed. The detailed form of the computation is shown 
in tahle XII of part III and differs from that used for 
single-stringer structures only insofar as necessary to take 
into account the fact that there is a sheet and a stringer 
on either side of the stringer boing adjusted instead of 
only one sheot ant st-rin^cr on one side. 
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Columns 1, 2, and 3 of tlMe XII 6 ive the values of 
the s triage* stresses , Oq , and cr-g • They are listed 

in t : .iis sequence so as to separate the values of &a and 
CJq , which remain constant during the adjustment of string- 
er 3, iron the stresses ojg and the other quantities 
that change during the adjustment. 

Columns 4 to 7 give the computation of the shear force 
in the panel between stringers A and 5; all properties 0: 
this panel are denoted by the superscript ' A3. 

Columns 8 to 11 give the computation of the shear 
force in the panel between the stringers 3 and C; all 
properties of this panel are denoted by the superscript 3C. 

Column 12 gives the difference between the shear 
forces in the two panels for each bay 

A3 3C 

D = AS - AS 

CS CE 

Column 13 gives the force = QkJL*« 

Column 14 gives the increments AF-g . 
Column 15 gives the shear fault 

S? = D ~ AF E 

Columns 16 to 19 give the shear-f ?Ail t correction stress 
A °3 in - r: " lo ee r with the columns 11 to 15 of the single- 
stringer* computation. 

After several cycles - say five cycles - the adjust- 
ment^ of stringer B is stopped, and the stringer is locked 
at the stresses thus obtained. 'The process of adjustment 
has again upset the static equilibrium; that is, the ex- 
ternal force at any cross section will not "be exactly bal- 
anced by the summation of the internal stringer forces as- 
sumed to exist at this stage. Static equilibrium is re- 
stored as before by increasing the stresses in the stringer 
that will be adjusted -next, nanel~, stringer C. 

Stringer 0 is now unlocked and adjusted, and the 
procedure of adjusting and restoring equilibrium is con- 
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tinuod until the center stringer is reached. The entire 
process is then repeated several bines until successive 
values of all stringer stresses in the structure ar^ in 
suf f i ciently close agreeneiit • 

In the case of a cambered cover, it is necessary to ' 
introduce the same modification as discussed for single- 
stringer beams, "based on the modified "basic equation (ic 1 ) • 
After column 4 of table XII a column must he added for 



Principle and scope of method.-' The mo st convenient 
and the most rapid method of analyzing structures with cut- 
outs is the indirect, or inverse, method. The analysis by 
the indirect method is made in two steps. First, the struc- 
ture is analyzed for the "basic condition that exists "before 
the cut-out is made. The results of this "basic analysis 
are used to calculate the internal forces that exist along 
the "boundary of the proposed cut ~ out . External forces equal 
and opposite to these internal forces are then introduced; 
these external forces reduce the stresses to zero along the 
"boundary of the proposed cut-out, and consequently the cut- 
out can now "be made without disturbing the stresses. 

The external forces introduced to reduce the stresses 
along the "boundary of the cut-out to zero will "be called 
the "liquidating" forces, a term used by R. V. Southwell in 
& somewhat different moaning . In general, it will be im- 
possible to calculate accurately the stresses that these 
liquidating forces set up at a distance from the cut-out. 
Some simplification of the problon is permissible, bccaf.se 
the liquidating forces form self-equilibrated systems so 
that, by St. Vcnant*s principl e , their effects become neg- 
ligible at some distance from tho cut-out. In order to ob- 
tain numerical answers, however, it is necessary to make 
very stringent simplifying assumptions, .and the method can 
therefore be applied only to reasonably small cut-outs. 

The treatment given here is confined to structures 
having distinct stringers* For cases in which the string- 
ers and the shin arc fused into a homogeneous unit, it is 
preferable to use the standard methods of the theory of 




similarly, after column 8 a 
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elasticity, and none solutions of the cut-out pro olen for 
such cases may be found in publications on the theory of 
olasti city. 

gf foots ©f refcOTing a ^kip |>a^el^ Figure 17(a) shows 
the internal shear forces that exist along the edges of a 
skin panel hounded by two stringers and two ribs. The di- 
rections of the force arrows are the positive directions 
in accordance with the general sigE conventions. In order 
to reduce the shear stresses along the edges of the panel 
to zero, external or liquidating shear forces are intro- 
duced as shown in figure 17(h), which are equal and oppo- 
site to the internal shear forces; only the forces acting 
on the main structure are shown in figure 17(h) "because the 
stresses in the skin panel itself are of no interest. 

In most practical cases $ the stringer areas and the 
skin thicknesses just outboard of the cut-out are the sane 
as those just inboard of the cut-out. The stress- 
distribution set tip by the liquidating forces will then he 
symmetrical about a chordwise line "bisecting the cut-out* 
Figure 17(c) shows schematically the stresses sot up in the 
stringers with the signs appropriate to the case where the 
basic stresses are positive. The figure indicates stresses 
only for the two stringers 'oorAcrin^ the cut-out; the 
stresses in the other stringers are snail enough (as will 
be shown experimentally in pt . Ill) to be neglected in view 
of the fact that the changes in stress distribution caused 
by a small cut-out arc snail compared with the basic 
stresses. 

The assumption that the liquidating forces of figure 
17(b) set up stresses only in stringers C and D is equiv- 
alent to assuming that the skin panels 30 and DE are 
rendered inoperative by slotting them lengthwise. Under 
this assumption, the problem becomes identical with the 
problem of the free panel shown in figure 18. The analyti- 
cal solution for the free panel is given in reference 1; 
for the present purpose it can be simplified by assuming 
that the structure is very lon 0 on either side of the cut- 
out. The forces in the stringers inboard and outboard of 
the cut-out are then given by the formula 

? = ± § T Q t * (18a) 

where T Q is the basic shear stress existing in the panel 
before the cut-out is made. t is the thickness of the 
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panel, L Is the lengt.h of the cut-out panel, and 
the sliear^lag parameter define d by 



Gt /l 
Ed \A C 



K 2 . * ±) (1ST)) 



The signs of the stringer stresses set up by the 
liquidating forces ? are indicated in figure 17(c) for 
the case of a positive b a sic shear stress t . 'The shear 

stresses set up hy the liquidating forces are given by 

T a I T 0 KLe (18c) 



and are of such a direction as to increase the basic shear 
stresses. Within the region of the cut-out, the stringer 
forces vary linearly between the maximum values obtained 
by setting x = 0 in formula (18a). The convention for 
measuring x in formulas (15a) and (13c) is Shown in figure 
17(a), 

The shear stresses given by formula (18c) are probably 
conservative because some of the shear load is taken by the 
adjoining panels, which are assumed tc be inoperative in 
this simplified theory. Conversely, allowance must be made 
for increased shear stresses in the adjoining panels. Con- 
st Aerations of continuity indicate that, in the immediate 
vicinity of the corners of the cut-out, the maximum shear 
stresses .in the adjoining panels BC and DS of figure 17 
should be taken as equal to the maximum stresses given by 
f o rmul a ( 1 8 c ) , 

jgC&ftfcg fff cu tting strinffeys»> Figure 19(a) shows a 
cut-out obtained by removing three skin panels and cutting 
two stringers. The effects of removin^; the skin panels can 
be calculated b f ^ the method described in the preceding sec- 
tion. The effects of cutting the stringers are represent- 
ed by the liquidating forces P shown in figure 19 (a) . 
The liquidating forces cause compressive stresses in the 
cut stringers and tensile reactions in the uncut stringers 
if the "basic stresses are positive, that is, tensile. By 
analogy with the preceding case of the skin panel, it mav be 
assumed that the tensile reaction to the liquidating forces 
is entirely furnished oy the two stringers bordering the 
cut-out; the stress system shown in figure 19(b) is cased 
on this assumption, and the numerical solution is obtained 
by considering one cut stringer and the adjacent continu- 
ous stringer to work together as a free panel. 
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Tho solution for tho froo panel (fig< l g ) of infinite 
length is 

^ a W Q * Pc~ Xx (19a) 

I cj, = p./a, a 2 = P 2 /A 2 (19b) 



with K defined by 



If Gyr;;:ctr" about a longitudinal lino through the' center 

of the cut-out is assuned, the nunorical solution for the 

cut-out is obtained in the first approximation setting 
in formulas (19 b) and (I9d) 



A 1 - Ajj = A-^ A 2 = Ap =a Aq. d = b (20) 

v7hore b denotes temporarily the off cctivc half-vridth of 
the cut-"-i-.t. The tests to oo described in part III indi- 
cate, hov;ever 5 that, even vs-boii only one stringer is cut, 
it is justifiable to assure that seyeral of the contirmouc 
stringers participate in furnishing the reaction to the 
1 i nui dating forces. The simplest assumption that can be 
nado about the participation of other stringers is ox- 
prcssod h^ setting 

A 2 . A & + A R o" i/% + A ie " 2d/l3 + . . . (21) 

uhor. formulas (19) are used. The stressor, caused "o' r the 
11 qui dating forces are then 

cr H = cr 2 o~ A /"° (?-2) 
o-j m cr 2 G -2d/'o 

Wb.cn onl-^ one stringer is interrupted, half of it is con- 
sidered as constituting A lt 7Th.cn n stringers are in- 
terrupted, the n/2 strii-i^ers on each side of the cut-out 
are considered to constitute A^, and the? are assumed to 
be concentrated at their cor.r.or. centroid to deterninc tr. 
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It is apparent that the use of f omnia (20) will "be 
conservative for stringers B and G and the skin panels 
oct',7con them but s brie what unconscrvat i vo for stringers and 
panels distant f ron the cut-out. 

At present, insufficient theoretical or experimental 
knowledge is available to define the limits within which 
the method presented here may he safely uccd. It would 
seem advisable to consider this method as giving only a 
first approximation when more than three stringers are in- 
terrupted oy the cut-out. The method of shear-faij.lt re- 
duction must he resorted to in such cases to improve the 
accuracy of the results. 

II. EXPBSI. aiSJ Ai VERIFICATION CP 
i a 1 THEORY 01 SHE A H LAG 
TEST OBJECTS AND TEST PROCEDURE 



NotT HA OA .tests . - Previous experimental investigations 

on shear lag have "been generally confined to panels and 
"beams of constant cross section; it was therefore consid- 
ered desirable to obtain experimental verification on a 
"bean \rith a variable cross section. Although the cross 
section can "be varied in a number of -.va~ r s, it was deemed 
most important and instructive to verify the influence of 
tapering; the cross-sectional area of the stringers. 

A skin-.st ri n^er panel was therefore "built as sho\m 
in figure 20 and tested in throe different set-ups. A pho- 
tograph of the second set-up is shown in figure 21 # In 
order to "obtain a sensitive check on the theory, the panel 
was designed for large shear-lag effects by using v a lar£e 
ratio of stringer are?, to sh«^et area. 

The tension panel was then converted into a beam by 
adding shear vebs; a cress section of the beam is sho\7ii in 
figure 2-8,1 and figure 2.? shovrs the inside of the beam with 
strain ^ages set up at one station. This beam is desig- 
nated beam 1. B ear. 1 was also tested with t.-o small cut- 
outs and two lar^e cut-outs located svmr.c t ri call Tr to the 
longitudinal axis. Figure 24 shov/s a st rain^gage set-up 
on the beam with the large cut-outs. 

After the cut-out tests rcro completed, the beam was 
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cut off just outboard of the first bulkhead, producing a 
very short, v:idc bead designated bean 2 . The test set-up 
for this short "bean is shown in figure 28 # 

It was also considered desirable to verify the valid- 
ity of the theory In the Uniting case of a hean without 
stiffeners. The dimensions of .a oean "ouilt for this purpose 
designated bean 3, are given in figure 2§, and the test 
set-up is shown in figure 27. In order to obtain a sensi- 
tive check on the theory, the hear. was r.ade quite short. 

As indicated in figure 26, hean 3 was tested in two 
conditions: first without corner flanges (original cross 
section) and then with corner flanges consisting of flat 
strips riveted to the cover as close to the corner as pos- 
sible (modified cross section). 

The hean was ouilt and loaded syrmc t ri cally about a 
transverse plane; it was thus possible to realize the con- 
dition of a built-in end and at the sane tine to neasurc 
strains directly at the root section. 

All strain readings were taken with 2-inch Tiickernan 
gages. These gages were always used in pairs on opposite 
sides of the sheet or the stringer to elininatc as far as 
possible the effects of local bending* Tenperatiire vari- 
ations during the tests were confined to 1° F , limiting the 
error in stresses to about 50 pounds per square inch. 

The load was applied in four equal steps in all of 
the cases except one, in which case three steps were used 
(bean 1, case 4). The stress readings plotted correspond 
to the highest test load used but were obtained by drawing 
the best-fitting straight lines through the load-strain plot; 
and correcting for zero shift -hen necessary. The fric- 
tion of the loading apparatus was measured several tines 
during the tests and was found to be 2 percent, unless 
otherwise .noted on the spanwise stress plots. Corrections 
have been applied for friction. 

Youn 0 f s r.oduli for the stringers were deternined fro- 
several specimens cut from the beans after the tests had 
been conpieted. For the sheet used to manufacture bean ^, 
the modulus was determined fron several test coupons cut 
fron the sane sheet fron which the bean was fabricated. 
The r.oduli obtained are noted on the drawings of the spec- 
i mens . 



In all these tests the Duckling stress of the sheet 
was never exceeded enough to cause an appreciable reduc- 
tion in the average shear modulus. In r.any tests there was 
no visioie "buckling at all. 

QL& tea t s re-anal ed^ Because the net hods of analysis 
proposed in this paper are relative!:' new , it seens desir- 
able to "buttress then with as nany experimental verifica- 
tions as possible. An effort was therefore na&s to secure 
all available test results and to analyze then bjr the pro- 
posed methods. It was found, however , that many published 
tests vrere -of doubtful value for furnishing quantitative 
checks because very thin sheet that buckled; at low loads 
had "been used in these tests; the effective shear modulus 
could not, therefore, he calculated with sufficient accu- 
racy for a quant itative . check, The tests considered usa- 
ble were a test on a compression panel r.ade by Thite and 
Antz (reference 4} and two beam tests r.ade by Schapitz (ref- 
erence 5), The oear. tested b^ r Schapitz is shown schemat- 
ically in figure 28. 



TEST RESULTS AITD COMPARISONS ~ITE THEORY 



i^ e thgds_of „analv_si s_used. All calculations vrere r.ade 
Tof analyzing the substitute single-stringer structure by 
means of the recurrence formula. The stresses in the 
stringers were computed by using the method of chordwise 
di st ri "oution as described in part I of this paper, includ- 
ing the correction for a finite number of stringers. Un- 
less otherwise noted the calculated rqsuits showg as curves 
in the figures are those obtained with the second approxi- 
mation for the substitute \ridth, 

Part I does not b ive explicit rules for determining 
the width b 3 of the idealized sheet "ootweoa stringers 
when the st. ringers are arranged as in hear. 1. The sheet 
tms assumed to he clar.ped betweon the opposing stringers 
with an effectiveness of 50 percent; in other words , the 
calculations were r.ade as though the stringers '.-ere attached 
"b"- two rows of rivets separated by half the width of the 
s t ringer s ♦ 

Mew H^QAr. test s.- The panel was tested under three con- 
ditions, as schematically indicated in figure 29. Figures 
30 to 32 show the experimental and the calculated results 
in the form of spanwise plots of stress. Figures 33 and 34 
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show the corresponding chordwise plots for the first two 
c a ses « 

The agreement between oxporir.er.t and theory is very 
Sati sf actery except near the root in cases I and 2. I&o 
c xporir.cntal points in this region scatter badly about a 
Dean line (figs. %Z and 34). Integration of the measured 
s t r e s s c s over the cross section b ivcs internal forces that 
agree dtMn about 5 percent with the external load, indi- 
cating that the strain measurements are fairly accurate Vat 
that there was sor.;e irregular "behavior of the structure. 
It v;as thought that this irregularity night oe caused by 
play in the bolt holes at the root; several holes were 
therefore carefully reaped out for the next larger size of 
holts before raking the hear: tests, and the chordwise plots 
of stresses for the beans uere much more regular. 

3ea:.: 1 was tested under the four loading conditions 
shown in figure 15. The spanwise stress plots are shown in 
figures 36 to 39 • The agreement between tests and theory 
is very satisfactory for the most highly stressed stringers 
near the flange and for the flanges themselves, except for 
the fact that the experimental stress in the flange at the 
station nearest the root is slightlv high in cases 1 and 4. 
In the stringers near the center line, the experimental 
stresses are higher than the calculated stresses near the 
root in cases 1, 3, and 4. It is "believed that the dis- 
crepancy can probably he charged to the assumption that 
the sheet was 100 percent effective in contributing to the 
stiff enor area. There are fairly consistent indications 
from a numher of tests that this assumption is too optimis- 
tic when the ratio Oy/cr L is large. A similar observation 
was made in reference 5 . This remark applies both to the 
compression side when the stresses are helow the buckling 
stresses for the sheet and to the tension side. On the com- 
pression side, the well-known effective width of the sheet 
nust oe use.d when the sheet has buckled. 

The results on bean 2 are shown in figure 40 # In view 
of the fact that this bead has .an extremely small ratio of 
length to width as well as a small shear-lag parameter K, 
the agreement is excellent. 

Figure 41 shows the results of test 1 on bear© 3. Be- 
cause the beam is s^nnetrical about the longitudinal axis as 
T7C11 as the transverse axis, there are four stress values 
for each station. It will be noted that in most cases the 
four values agree ver- closely, which indicates that the beam 
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slicked excellent • cyr.notr;; of strain ahout both axes. 

This feast is a rather crucial tost or. t ho range of 
validity of the theory. It has "boon held o ,r sor.o investi- 
gators that the theory of shear lag as developed in this 
paper would not ap^ly to the Uniting case where the ele- 
ments of the co Ycr-carr"ing shear (the shoot panels) and 
the elements carrying normal stresses (the stringers) are 
no&ged into a single unit, namely, a shoot. Figure 41 
sho-s that this opinion is too pessimistic; the a G rccnoi:t 
is not perfect, out the maximum flange stresses, which are 
of paramount interest for design t are predicted fairly well. 

The main difficulty in applying the theorv to the 
case just discussed lies in the fact that Ay "becomes very 
sr.all compared with At; the flange area consists only of 
the area l/o ht , which expresses the participation of the 
shear web in the bending action. For snail ratios of 
to Al, the chear-lag parameter E becomes ver^ lar^e 
and sensitive to errors in Ap . The difficulty is ooviat- 
ed when a corner flange of reasonable area is provided; in 
"built-up structures, such a cornor flange is usually pro- 
vided in the form of an angle for riveting the cover to 
the shear -eo. In "beam 3, a corner flange was provided by 
riveting flat strips aIon b the ed^es, as shown on the sec- 
ond cross section in figure 25. The test results for this 
condition are plotted in figure 42 and show excellent agree- 
ment with the theory. 

01 d_ t _e s t_s . - Figure 43 shows the experimental and the 
calculated results for the compression panel described in 
reference 4. 

Figure 44 shows the results of the test on the "beam 
described in reference 5 for a load applied &t the tip. 
Figure 45 shows the tost results for the s am o bean under 
loads distributed as indicated in figure 28. The agreement 
is fairly satisfactory. 

Gut -out tests.- The approximate method of analyzing 
cut-outs described in this pap^r is based on the assump- 
tion that a pair of equal and opposite forces applied to 
adjacent strin 0 ers does not affect other stringers very 
much. A special test was made on beam 1 to verify directly 
the validity of this assumption. Two equal and opposite 
forces of ? - 1162 potl&&s were applied to bolts at the 
intersections of the rib at midspan with stringers D and E # 
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Figure 46 shows the experimental stresses and the stresses 
calculated under the- assumption that only stringers D and 
E are stressed. 

Figure 4.7 shows the results of the test on "hear. 1 
with snail cut-outs located as shovm *by the sketch; only 
the shin was cut out in this case. 

Figure 4-8 shov?s the results of the test on bean 1 
with lar^c cut-outs located as shown "by the sketch and in 
figure 24. The a^reer.ent between theory and experiments 
for the cut-out tests is very satisfactory except for the 
discrepancies already noted in the tests on the sa.no bean 
without cut-out s . 



III. HUMERI CAL EXAMPLE S 
IDEALIZATION OF CROSS SECTIONS 



Profrlgn To find the idealized cross section of a 

bean with open-section stiff oners: 

The actual cross section of the bean is shown in fig- 
ure 49(a). The effective width of the sheet for nornal 
stresses is to be' taken as w = 20t . 

The idealized width & for shear deformation (fig. 
11) is equal to the width between rivet rows, that is, 4 
i nche s . 

The area of the idealized flange is obtained by add*- 
Ing the following areas: 

( sq in . ) 

Corner afcgle 0,300 

Skin fron corner to rivet line 

(0.375 X 0.040) ... . , 015 

Equivalent of web (l/o x S.*00 x 0.065 ) . . .065 
Area of idealized flange 0.380 

The first stringer innediatelv adjacent to the flange con- 
sists of onl- the effective width of skin; the area is 

20 x 0.040 x 0.040 = 0.032 square inch 
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Each of tie next tvo stringers consists of a stiff encr and 
a douole strip of skin; the area of each idealized stringer 
i s therefore 

A f 0,200 + ?. x c?0 x 0.01-0 x 0.040 = 0.254 square inch 

The stringer at the center line has one-half of this area, 
or 0 # 132 square inch. 

The total area of the longitudinals is 

A L $ 0.032 -f 0.264 + 0.254 + 0.132 = 0.692 square inch 

The idealized cross Section is shown in figure 49(h). 

ProMem 2 • - To find the idealized cross section of a 
"bean v:ith closed-section stiffeners: 

The actual cross section of the oearri is shovm in fig- 
ure 49(c). The effective ri&th of the sheet is to be 
taken as w = 20 t . 

The effective width o 2 for shear deformation is, oy 
formula (18 ), 

o = iiS- = 0.75 inch 

6.040 x 3,00 

The idealized width from the flange to the first stringer 
i s therefore 

d 5= 3.2 5 ♦ I X 0.75 = 3.53 inches 

and the idealized 'vidth of the second and third panel is 

d = 2.50 + 0.75 ~ 3. 1 25 inches 

The areas of the flange A ? and of the ^"irst snail string- 
er are the sa::e as in /pro olen 1. 

The area of the second as well as of the third ideal- 
ized stringer is obtained "by adding the following areas 
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Csq in, ) 

Hat section . . 0.250 

Skin ootwoon rivets (1.5 x 0.040) .050 

I wo strips of skin (2 x 20 x 0.040 X 0.040) .054 

Area, of idealized stringer ; 0.384 

The stringer at the center line has one-half of this area, 
or 0.19 2 square inch. 

The total area of the longitudinals is 

A L = 0.032 + 0.^84 + 0.384 + 0.192 = 0.992 square inch 

The idealized cross section is shown in figure 49(d). 

ANALYSIS 01 A IvlULTISTRIKQER BEAH f OBTAINED 3Y THE USE OE 
THE SUBSTITUTE S I1TGLE-STRIHGER METHOD AND 
THE RECURRENCE FOSMTJXA 



Griven &ata.- Ei^ure 50 shows the idealized forn of a 
"bean; the problem is to find the stresses in this he am un- 
der the load* indi cat ed "by the use of the substitute single- 
stringer method and the recurrence formula. 

This idealized "bean is ver3 r nearly identical with the 
idealized forr. of "beam 1 discussed in part II. The follow- 
ing simplifications have "been nade : The slightly -tapering 
effective width of "bear. 1 has "been replaced "by a constant 
width; the slightly tapering effective depth, with a dis- 
continuity at the nidspan, lias, ooen replaced "by a constant 
depth; the load has "boon located exactly at the tip instead 
of at the actiial location of 0.55 inch iron the tip. ITone 
of these deviations ar.our.ts to r.oro than 2 percent- at any 
point, so that the results obtained in these numerical ex- 
amples can be compared quite closely with the correspond- 
ing calculated curves shewn in part II. 

From the data given in figure 50 , table I has "boon 
prepared to b ivo the data in the form required for the anal- 
y sis. 
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first ammrnxig&t 1qs^M._^ stJAait. .siSSiSrSiiiSfiSJe 
.structures*- The first appro xir.at i on to the substitute 

single-stringer structure is obtained "by combining the 
.stringers constituting A L . into a single stringer located 

at the centroid of A L . As indicated in figure 50, this 
•centroid is located 6,28 inches from the flange, and this 

distance is by definition the substitute width in the first 

approximat i on. 

The computation of the coefficients required for the 
analysis of the substitute beam is shown in table II. The 
values of Ay and A L are the same as for the actual 
structure and are obtained from table I. The shear-lag 
parameter K is calculated by formula (4). The substitute 
width bg just found is used' where b appears in this 
formula, so that 

Eos 5 * 23 

The coefficients p, q, and V are calculated by formu- 
las (3a), (3b), and (3c); because 6 and t are constant 
in this particular beam, the common factor G-t has been 
omitted from all coefficients. 

TTith the coefficients computed in table II, the sys- 
tem of equations for the X-forces (first approximation) is 
written in conformance with equations (5). The boundary 
conditions arc X 0 = 0 and V r + X = 0, 

- Xx (0.1400 + 0.1383) + X 3 (0. 1132) = -65 . 7 + 55 .5 

X^(0^1182) - X 2 (0. 1388 -1-0.1376) + X 3 (0. 1130) = -55.5 +65,3 

X 2 (0.1190) - X 3 (0.137S + 0.1352) +X 4 (0. 1191) =-55.3 +65.1 

X 3 (0.119l) - X 4 (0. 1362 + 0.1358) + X 5 ( 0 . 1200) * -55.1+ S5.0 

X 4 (0.1200) - X B (0. 1353 + 0.1347) + X 6 (0. 1201) * -55.0+65.0 

X 5 (0.120l) - X 6 (0.1347) = -55.0 

Those equations are then solved, and table III £>ives the 
final computation of the stresses in the substitute beam as 
obtained by formulas (9a) and (9b). 
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fteso&ft a^roxin&frtoa to smji ata*»t_ g^M 
atructu,re»«« The calculation of the second approximation 
■begins with the last two coluwis of laWe III. The param- 
eter Yb is obtained from figure 13 for each station, and 
the average value of Yb is computed. From figure 15, 
the value of 1 - (^v/la) corresponding to this average 
value of Yb is read, and the second approximation to the 
substitute width is obtained bv formula (17a). Actually 
it is no t necessary to compute the second approximation of 
bg ; it is possible to proceed directly to the new values 

of the shear-lag parameter K by dividing the values of K 

given in table II by the expression J%>{X ~ C^t/^ )j • 5Nt* 

ble IV gives the new values of X and the computation of 
the new set of coefficients p and q for the recurrence 
formula. lable V gives the values of the X-forces and the 
final stresses in the beam for the second approximation. 
As a check, the average value of Yb is again computed, 
and the corresponding value of 1 - (y^/h) is found. The 

factor J 2 [l - (r^/b)} differs by onir 1 percent from 
the factor obtained in the first approximation; the second 
approximation may therefore be considered the final a,pproxi 
mat i on . 



0^1culati^_of Ol^ordwise dirtrl'bution oS stresses,- 
After the final approximation to the stresses in the sub- 
stitute beam has been computed, the chordwise distribution 
of the stresses in the actual beam can be found. As an 
example, the calculation will be shown in deto.il for sta- 
tion 5 . 

According to table V, (7-, = 5000 pounds per square 
inch and 0*^ = 2573 pounds per square inch for station 5; 
the ratio ^Z 0 ! = °-535, and the corresponding value of 
Yb = 1.750 from figure 13. This value of Yb is entered 
in table 71, and the values of Yy for the two intermedi- 
ate stringers B and C are calculated "bjr proportion and 
entered in column 2 # - Next, the hyperbolic cosines are en- 
tered in column 3. The stress in the center stringer D 
is nor; calculated bv formula (16). 

5000 

CTqt = = 1573 pounds per square inch 



^.nd entered in colur.n 4. The stresses in the stringers B 



&nd 0 arc then calculated by f*o rt&ula (14) and entered in 
column 4. 



Column 5 ^ives the cm s s- so ct i on-rl areas of the spring- 
ers A S 4.. and column 6 gives the internal forces a A^.- 
The QVMi of these forces wi 11 not equ.al the force CT^i^. on 
account of the finite number of stringers, and a correc- 
tion nust be applied to all of the stresses 0" except to- 
the stress in strin 0 or A; the stress in stringer A must 
necessarily remain equal to • 

The correction is made as follows* The force F* is 
c L A L a 2673 x 0.771 s 2050 pounds. The force in stringer 
A is 140 pounds, as shown in column 6; the total force 
that must he applied b?' the center stringer D and the two 
intermediate stringers 3 and 0 is therefore 2050 - 140 = 
1920 pounds. The summation of the internal forces in the 
three stringers 5, C, and D as given in column 5 is only 
1715 pounds; the stresses 0" given in column 4 must there- 
fore he multiplied "by the factor 1920/1715 = 1.120 to 
obtain the final values of the stresses cr, which arc 
listed in column 7. As a check, the internal forces are 
again computed with the corrected values of cr; the summa- 
tion chochs exactly with the force Ft a 2050 pounds. 

The calculation of the chord\risc distribution of 
stresse.s is made in the same manner for each station; the 
results of the calculations arc given in table VII. 

ANALYSIS 0? I-IULTI STRINGER BEAM WITH CUT-OUT 



It will "be assumed for the example of a mult i st ringer 
"beam -with cut-out that a cut-out is made in the "beam shorn 
in figure 50 and analyzed in the preceding example; the 
skin panels A3 and 30 and the stringer 3. are removed 
"between stations 3 and 4^, corresponding to' the lar L :e cut- 
out in beam 1 described in part II. The effects of making 
this cut-out are to oc found. 

Iff ect s JPf tpmrim gfetn , qanejL A3.- The total shear 

force in the shin panel A3 between stations 3 and 4$ is 
found "by statics with the stresses given in table VII; it 
is equal to the s\im of the forces in stringers B, C, and I) 
at station 4-£ minus the sum of the forces in the same 
stringers at station 3. The result ef the simple calcula- 
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t i on is T Q tL = 494 pounds. The next stop is the calcula- 
tion of the parameter K "by formula (18b) • In this ease , 
the cut panel is "bounded by stringers A and B ; the areas 
Aq and A^ of formula (18 o) are therefore replaced 07 
A^ and A-g . In order to "be consistent with the assumption 

that the structure is the sane at the two ends of the cut- 
out, tho values of A^ rnd A3 used will "oe those valid 
for the middle of the cut-out. Formula (13b) gives there- 
f o v e 

= M^ISMZ (-1^ * *XJ) m 0.01078 
3.625 Vo.355 0*271./ 

f » 0.1038 

After these preliminary calculations, the solution can "be 
carried out in tabular form as shown in ta"blo VIII.. The 
value of ? at stations 3 and 4-§ is f T Q tL = 242 pounds; 

at tho other stations P - 242 o"^ :z pounds according to for- 
mula (lSa). The calculation cf the stresses P/A^ and 
P/A-g i s self -explanat o ry . 

Iffecta of ,r^Q7i^' skin .BO;* For panel 30, 

the shear force is found by subtracting the internal forces 
in stringers C and B at station 3 from the forces at sta- 
tion the result is 

T Q tL be 19 5 pounds 

The value of K is found from 
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= 0.40_x 0^015 + „.1„ X ) * 0.012; 

3.625 \0.271 0.271, 



I = 0.1104 

Table VIII shows the details of computing the stresses 
P/A 3 and P/Aq caused by rerovin^: the skin panel 30. The 
last four rows of the table give the stringer stresses in 
the beam, obtained 0^ superposing on the stresses of table 
VII "the stresses caused tfrjr removing the two skin panels A3 
and 30. The si b ns of the stresses arc determined ov com- 
parison with figure 17(c); at station 3, for instance, 

0" B 3 2370 + 945 - 381 = 2S34 pounds per square inch 



44 



where 945 pounds par s guar. 6 inch is the stress caused "by 
renoyi ng panel AB , and 381 pounds per square inch is the 
stress caused "by renoving panel BO. 

If feet q£ cutting ,gt ringer.- According to the stresses 
listed in table 71 II, the stress in stringer 3 at station 
3 is Og =2934 pounds per square inch. The internal force 
at the outboard end of the cut-out is therefore 
2934 x 0.256 s 752 pounds. At the inboard end of the cut- 
out, the force is 2514 x 0.287 = 750 pounds. The region 
around the cut-out is nov; divided into four free panels so 
that formulas (19) can be used. Two of these panels are 
inboard of the cut-out; for the first panel 

A X = I A 3 a - d A 2 = A A 

for the second panel 

&l - i A 3 and A 2 = Aq + e A s • 

by z o rnul a ( 21 ) , all areas bein^ those at station 4§. For 
simplicity, it will be assured that the tvro panels have 
the sane shear-lag parameter II, and E will be computed 
by usinj; the average of th- tv;o ^iven values of A 2 . The 
result is 

** * &M~lL£*m f-JU, + -J. ^ = 0.01518 

3.625 \0.35 3 0.1435/ 

E = 0.127 

for tlxo inboard panels. 

The other two free panels r-re outboard of the cut-out 
and are defined in the sane r.ianner; the calculations are 
made vrith the areas at station 3. The shear-la£ parameter 
i s (:± von by 

r * = OjM^^im (L-l_ + -JL> = 0.01813 
5.625 VD.71S 0.128/ 

X m 0.1547 

The cnlculnti on itself is given in table IX. The stresses 
caused Igfc cutting otyiagfif* 3 shovrn in this t-olo "re 
superposed on the final stresses shown in tnhlc VIII to 
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Obtain the final stresses in the stringers. The stresses 
in stringer D caused "fay cutting strin b cr B aire obtained 
by formula (22) as e~ 1 P/A 2 . 

When the results of this computation are compared with 
the curves in figure 43, it should he borne in mind that an 
additional snail correction must be made for the actual 
test "because removal of the skin panels reduces the areas 
A and A^ in the region of the cut-out. 

ANALYSIS 3Y SUCCESSIVE SHEAR-FAULT REDUCTION 



Anal y si a of single- st rin&e r beam . - The method of ana- 
lyzing a sin^le-st ringer beam by successive shear-fault 
reduction will he demonstrated on the substitute single- 
stringer hear, analyzed previously by the recurrence formu- 
la. The basic data for the beam are those ^iven in table- 
I; for the substitute width, the second approximation 
bg = 5 .28 x 1.090 = 5.85 inches was used. As the initial 
assumption, the stresses in the flange were arbitrarily 
assumed to be 1.40 tines the stresses given by the Mc/l 
formula. Table X gi *e s the first cycle of the computation; 
a comment on the form of the computations is ^iven in part 
I of this paper. Table XI gives the second cycle of the 
computation, starting with the values of cr ? found at the 

end of the first cycle . As a general check on the compu- 
tations, the sun of the shear faults is 'shown for both cy- 
cles! it \rill he noted that it has decreased from 843 to 
764 pounds. 

Analysis ..of . nultistringor As afc example for the 

analysis of a mult i st ringer "ream, the bean of fi^re 50 is 
again used, and a typical cycle of adjustment for stringer 
3 is shown in table XII. Because the example is illus- 
trative, the stress values <X^ , CTg, and 0"^ \ierc not 
assumed arbitrarily but were taken from table VI I 9 the fi- 
nal result of the previous analysis. The shear faults are 
therefore very small, and the adjusted stresses cr^ o.rc 

practically identical vrith the initial stresses. The small 
differences that exist arise from tvro reasons. The first 
reason is the' limited numerical accuracy of the process. 
This numerical accuracy is determined b: r t\?o factors: the 
number of bays used and the accuracy of mul t ipli cat i on and 
division. These operations vf 8 r e carried out >ri t h a 10-inch 
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slide rule in all numerical examples given in thin report. 
The second reason for the failure of table XII to show ex- 
act agreement between the initial and the final values of 
(T-p lies in the slight differences between the basic as- 

sumptions. The recurrence formula is "based on the assump- 
tion that the cross section is. constant in each hay, out 
the stresses var^ nonlinearly in each hay* The shear-fault 
reduction method, on the other hand, assumes that all 
stresses vary linearly in each oay. 



CONCLUSION 



The theory of shear-la> action presented in this paper 
is "based on the concept of idealized structures consisting 
of stringers carrying longitudinal stresses, of sheet carry- 
ing shear stresses, and of transverse ribs .infinitely close- 
ly spaced and of infinite stiffness. The test results in- 
dicate that this theory is acceptable as a oasis for prac- 
tical stress analysis , because, in b eneral, the. differences 
between test results and calculated results in the criti- 
cal regions are smaller than occasional scatter of test 
results caused by uncontrollable irregularities in tho be- 
havior of the structure. 



Lan b le* r Memorial Aeronautical Laboratory, 

national Advisory Conr.it tee for Aeronautics, 
Lan 0 le^ Field, 7a. 
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APPENDIX A 
SYMBOLS 

A ci oss-section?.! area, sq in, 
I Young;* s modulus , It)/ sq in. 
F internal force, 11 

G effective shear ao dulus * lb/sq in. 

I geometric moment of inertia, in. 4 

M shear-lag parameter (equation (4)) 

L length | in. 

M belt ding moment, in. -Id 

P external load, lb 

S shear force, lb 

SF shear fault (equation (SS-3)} 

SFC shear-fault correction (equation (SS-4)) 

Y auxiliary parameter (equation (14)) 

1 half-width of structure, in.; with nunerical sub- 
scripts, distance between stringers (fig* 12), in 



b' developed width, in. 

h depth of bean, in. 

t t hi cknes s , i n 0 

w effective width 

x distance parallel to center line 

y distance fror. center line 

» 

y shear strain 
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0* direct (normal) stress, lb/sn in. 
T shear stress, lb/ s q in. 

T Q basic shear stress existing before a cut-out is nade , 
1 b / s q i n « 

Superscripts have the following significance 8 

P theoretical values based on the assumption that plane 
cross sections remain plane 

Subscripts have the following significance: 

C cover sheet 

32 external (applied) 

J flange 

L longitudinal 

S substitute 

st stringer 

T total 

W 'shear v:eb 

CE occurring in the cover sheet and obtained by the 
elastic relation. 

CL center line 

i in "board 

o outboard 

av average 

e effective 



49 



APPENDIX 3 
CC:.:?ARIS0:7 B10PWSEK DI7PSR31T7 SOLUTIONS 
0? THE SHEAR- LAG PROBLEM 



The oasic shear— lag proolon is the pro Die::: of a box 
beam with constant cross section. In 1930 Younger pub- 
lished a solution of this problem (reference ?). In 1937 
the author published a slightly different solution, the 
constant-stress solution (reference l). In 1938 Reissncr 
published a third solution (reference 3). If the flange 
efficiency 7) of a box be a..: is defined b~ the ratio of 
the M g / 1 stress to the actual flange stress, all three 
solutions can be reduced to the sane form, nar.clv, 



tanh_P 



where F is a function of 
cal properties of the box. 
as folio w s : 

I r 1.371 b/L 



the geoneiri cal 
This function 



and the physi- 
7 is defined 



(Younger , reference 7) 



( Euha , refer e r. c e 1 ) 



Lei ssner , re: or once 



3) 



It qri-ll be seen that the three solutions are identical in 
forn and differ only slightly in the numerical constant. 

All three solutions involve seoe simplifying assunp- 
tions, and any one of the three could be used equall Tr well 
as a basis for building up approxinat e solutions for beans 
of variable cross section. All three solutions, however, 
lead to the result that the flange efficiency is constant 
aleng the span* A glance at figures 41 and 42 indicates 
that this result cannot be :.:ore than a rough approximation; 
the flange stresses on these figures are not straight lines, 
Tor this reason, the treatnent of the bean -vith variable 
cross section a s t>t c s cnt c d in this paper was not based on 
any of these solutions. 
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Of the throe basib solutions given, only Reissner 1 s 
solution is of Such a nature that the underlying a s sump-* 
t ions car. "be physically realized without difficulty (con- 
stant cross section, concent rat ed load at tip). At the 
tine of publication, it was stated that the solution is 
applicable only when the cover consists of corrugated 
sheet (reference 8); it was stated l£.ter (reference 9) 
that the solution applies also when the cover consists of 
a flat sheet. Reissr,er ? s solution is therefore shown in 
figure 41; it will he seen that, at sonc distance from the 
root, it is a fair approximation, out at the root the ex- 
perimental shear-lag effect is nearly twice as lar 6 o as that 
predicted by Rei ssner ! s solution. 

The series solution given by TTinn^ (reference 10) is 
"based on the sane principles as the solutions listed and 
is therefore open to the sane o ejection in that it cannot 
give norO than a very rough approxiraat i on . In view of 
this fact, the labor of using a solution by series is hard- 
ly justifiable. 

The solution given by G-oodey (reference 11) is iden- 
tical with the solution of the single-stringer "bean given 
in reference 1. G-o^de- also jives one case not included 
in reference 1, nanoly, the case of uniforr.lv distributed 
loading. 

A very complete and elaborate net hod of shear-lag 
analysis has been presented by Eonor and Koller (refer- 
ence 1?.). The idealized structure consists of stringers, 
sheet, and transverse ribs. The transverse ribs are fi- 
nite in nunber and of finite stiffness; the r.ethod is 
therefore nore complete than the methods presented in this 
paper. Cor.parat i ve calculations nade in reference 12, 
however, show that the rib stiffness has only a snail in- 
fluence on the stringer stresses so that the sir. 1 .?! if ying 
assumption of infinite nunber and stiffness of the ribs 
results onl T - in ver ir snail errors. This conclusion drawn 
by Ebner and Koller fror. their theory is anply confirneel 
by the good agreement between* "the experiments and the anal- 
yses presented in this paper. 

The nethod of reference 12 is rather difficult to fol- 
low; comparisons have therefore been confined to the anal- 
ysis of numerical examples given therein o rr the methods 
presented in this paper. The dimensions of the- structure 
analyzed in reference 12 are given in figure 51. 
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Figure 52 shows the results for load case 1, irhi c& 
is the case of an axially loaded panel. It will "be noted 
that in reference 12 there is given an "exact" method as 
well as an approximate one, the approximate net hod "being 
r econnended for. practice "because the exact method is quite 
cumber :;or.o . The solution nadc by the substitute single- 
stringer method agrees with the exact r.cthod of rcf er^ncc 
12 at all of the stations except one to the accuracy of 
reading the values from a snail graph. The r.axinun dif- 
ference "between the exact method and the present single- 
stringer nethod is only slightly larger than the differ- 
ence "between the two methods of reference 12 and is unim- 
portant for design purposes. 

figure 5? shows the results for the beam. The agree- 
ment between the solution of reference 12 and the single- 
stringer solution of this paper is very close except at 
the root,. where there is a difference of 3 percent on the 
flange stress and a difference of 20 percent on the stress 
in the center stringer* The agreement "between the solu- 
tion of reference 12 and the solution "by successive shear- 
fault reduction is good. 

It should he pointed out t-iat this numerical example 
represents the most severe test that. can possiol; r "be made 
of the powers of the substitute single- s t ringer method. 
The chordwise di st rioution method, which is an integral 
part of this method, is based on the assumption that there 
are infinitely manv stringers; the half structure analyzed 
here has only two stringers, which is not a very close ap- 
proximation to infinitely many stringers. 

The example may serve as a yarning, therefore, that 
in such extreme cases, the method of shear-fault reduc- 
tion should be used to refine the approximation obtained 
by the single— ringer method. From a practical point of 
view, the discrepancy found here between the method of 
reference 12 and the substitute single-rst ringer method is 
o.f little .interest because structures .with only two string- 
ers are not likely to be encountered in practice. 
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Tables l r 4 



jj}/E = 0.1+0; 
t = 0.015 in 



TABLE I 

BASIC DATA FOR ANALYSIS OF BEAM 

h = 6.18 in.; b = IO.875 in.; 





X 


A p 


A L 




II 


II 


-1 


Station 










n 






(in.) 


( sq in. ) 


(sq in.) 


(sq in.)' 


(Ib-in. ) 


(lb) 


(lb/sq in.) 


0 


0 


0.229 


0.502 


0.751 


0 


0 


0 


1 


8 


.256 


.556 


.812 


4,800 


777 


956 


2 


16 


.285 


.610 


.893 


9,600 


1554 


1740 


3 


24 


.510 


.665 


.973 


14, 400 


2331 


2395 


4 ' 


32 


•557 


.717 


1.054 


19,200 


3108 


2947 


5 


4o 


• 564 


.771 


1.135 


24,000 


3885 


3i^o 


6 


1*8 


.392 


.825 


1.217 


28,800 


4662 


3835 



TABLE XT 

Coefficients for Recurrence Formula (Second Approximation) 



Bay 


K 


KL 


tanh KL 


sinh KL 


P 




1 


0.0826 


0.661 


O.579 


O.710 


0.11^5 


0.1165 


2 


♦ 0786 


.629 


•557 


.671 


.1410 


.1172 


3 


.0750 


.600 


• 537 


.657 


• 1596 


.1176 


. 4 


.0718 


• 574 


.518 


.606 


.1556 


.1184 


5 


.0690 


.552 


.502 


.580 


.1374 


.1190 


6 


.0 665 


.532 


.487 


• 557 


.1366 


.1195 
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Tables 2,3 



COMPUTATION OF COEFFICIENTS FOR RECURRENCE FORMULA (FIRST APPROXIMATION ) 
fit 2 = O.OOO956 [j^ + jj-)- L = 8.00 taj TABLE II 



Bay 


Csq in) 


A L ^ 
(sq in.) 


1 

A p 


1 

A L 


1 ♦ 1 


K 2 


K 


KL 


ta nh KL 


Sinn KL 






At 


Y 


1 


0.21+2 


0.529 


1+.13 


1.89 


6.02 


3.00575 


0.0758 


O.607 




0.6I+5 


0.li+0( 


0.1177 


0.687 




2 


.269 


.583 


3.72 


1.72 


5-W+ 


.OO520 


.0721 


.577 


.520 


.610 


.1588 


.11°2 


.685 


66.5 


3 


.296 


.636 


3. 33 


1.57 


h»95 


.OQl+73 


.0688 


.550 


.500 


• 578 


.1376 


.1190 


.6P3 


66.5 


if 


.323 


.690 


3.10 


1,45 


u.55 


.001+35 


.0659 


.528 


.m 


• 553 


.1362 


.1191 


.683 


60.I 


5 


.350 


.7kk 


2.86 


i.5k 


14-20 


.ooi+02 


.0633 


.506 


.1+67 


.528 


.1353 


.1^00 




r b.O 


6 


.378 


.798 


2.65 


1.25 


5.90 


.00573 


.0610 


.1+38 


.453 


.508 


.131+7 


.1201 


. f> c >0 





STRESSES 
(1 



orruuoo 

[jb = 



IN SUBSTITUTE SINGLE-STRINGER BEAM 
y 

- r^) = 0.1+20 from fig. 15- 

D 



V 0.81+0 



= 1.090^ 



-for 



(FIRST APPROXIMATION) 



TABLE III 



Station 



(lb/sq in.) 



X 
(lb) 



(lb/sq in.) 



(lb/sq in. ) 



(lb/sq in.) 



(lb/sq in.) 



<*L 
OF 



Yb 



956 
171*0 

2395 
291+7 
51*20 

3835 



1+8 
113 
212 

373 
633 
105U 



188 
i+00 
681+ 
1107 
171+0 
2690 



111+4 
21I+0 
3079 
I+05I+ 
5160 
6525 



86 
186 
320 
520 
822 
1278 



870 
1551+ 
2075 
21+27 
2598 
2557 



D.760 

725 

675 

60C 

501+ 

392 
Total 

Average 



1.00 
1.10 

1.25 
1.50 
1.90 

2^0 
9.25 

1.5U 



Table JV_ 
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TABLE T 

STRESSES IN SUBSTITUTE SINGLE- STRINGER BEAM (SECOND APPROXIMATION ) 
from fig. 15. 

1 



Tables 5,6 



f(i . gt ) * 0.429 



y 2 (i . ik) 

D 



VO.858 



= 1.079 



3 



Station 


o p 

hVj/fln -fn } 


X 


X 

Ap 

flb/«a In 1 


Clb/«a in. ) 


X 

\ 

(lb/sa in ) 


(lb/sq in. ) 


OL 
CTp 


Yb 




1 


956 


ho 


156 


1112 


72 


88I4. 




0.795 


O.90 




2 


171+0 


96 


339 


2079 


157 


1583. 




.761 


1.00 




3 


2395 


187 


603 


2998 


282 


2113 




.70U 


1.16 




1+ 


2947 


355 


994 


3941 


U68 


214.79 




.629 


l.ljO 






3420 


576 


1580 


5000 


7^7 


2673 




• 535 


1.76 




6 


3835 


987 


2520 


6355 


H95 


261^0 


Tc 


.416 
tal 

rerage 


2.4c 
3.62 

Ml* 





TABLE VI 

COMPUTATION OP CHORDWISE DISTRIBUTION 
OF STRESSES AT STATION 5 



Stringer 


y 




00 sh Yy 


0 

(lb/aq in) 


A st 

(sq in. 


(lb) 


a 

llb/sc 4t) 


0 A g . 

(lb) 


£ D 


0 


0 


1.000 


1673 


0.149 




1873 


279 


C 




.587 


1.177 


1968 


.297 




2203 


655 


B 


§» 


1.174 


1.772 


2965 


.297 


881 


3320 


986 


A 




1.760 


2.992 


5000 


.028 


(i4Q) 


5000 


140 












.771 p.715 




2060 



U A = 2673 x O.77I = 2060 Correction factor = = 1.120i 



1920 lb 



Uncorrected values. 
Corrected values. 



TA3L3 VII 



Chordwise Distribution of Stresses la 3eaa 



Stat ion 


°A 

(lb/sq in. ) 


°3 ■ 
(lb/ so in. ) 


°C 

(lb/sq In. ) 


(lb/ sq in. ) 


1 


1112 


952 


840 


800 


2 


20 79 


1774 


1518 


1440 


3 


2S9S 


2370 


1936 


1790 


4 


3941 


2875 


2200 


1985 


5 


5000 


3320 


2203 


1873 


6 


6355 


339 9 


2020 


1530 



TABLE VIII 



EFFECT OF REMOVING SKIN PANELS 



Station 

X 


6 
12 


5 


4 


3 
0 


2 
8 


1 
16 


Effect of removing panel AB 


Kx 

p (lb) 
A^(sq lr 

p/a a 

( sq in 


.287 
69.5 

J .422 

165 

J .318 

219 


0.1*15 

.61+7 
156.7 
.392 

1+00 

.297 

527 


0 

1.000 
21+?. 0 
.377 

61+2 

.286 

81+7 


0 

1.000 
2I+2.0 
.331+ 

725 

.256 

91+5 


0.830 

.1+36 
105.6 
.305 

31+6 

^235 

1+50 


1.660 

.190 
1+6.0 
.276 

167 

.215 

211+ 


Effect of removing panel BC 


Kx 

e -Kx 

P(lb) 

P/Ac 


1.325 

.266 
25.9 
81 

81 


0.1+1+2 
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Figure 3.- Convention for coordinate axee . 
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Figure 11.- Idealisation of multistringer cross section. 




Figure 13.- »t*niar& ayafcols for width of panels. 




Figure 13.- Auxiliary graph for determining cliordwise distribution of stresses. 
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figure 19.- Effects of cutting etriDgers. 
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Figure 32.- Cross section of beam 1. Oorer of bona is panel 
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Figure 39.- Diagram of leading oaass for panel, F= 1300 pounds. 




Figure 24,~ Test set-up for beam 1 with cut-outs 
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Figure 30.- Stresses in panel, case 1. • 
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Figure 32.- Stresses in panel, case 3. 
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Figure 36.- Comparisons between calculated and experimental stresses in 
beam 1, case 1. 
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Figure 34.- Chordwise distribution of stresses in 'panel, case 2. 
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Figure 36.- Diagram of loading cases for bean 1. P * 600 pounds 

on seen shear web fox cases 1,2, and 3; P « 3*5 
pounds on saon shear wob for case 4. 
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Figs. 37,38 
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Figure 37.- Comparisons between calculated and experimental stresses in 
beam 1, case 2. 
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Figure 38.- Comparisons between calculated and experimental stresses in 
beam 1, case 3. 
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Figure 39.- Comparisons between calculated and experimental stresses in 
beam 1, case 4. 
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Figure 46.- Comparisons between calculated and experimental stresses for 
special tests on beam 1. 
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Figure 41.- Comparisons between calculated Figure 42.- Comparisons between 

and experimental stresses calculated and experimental 

in beam 3, case 1. Friction less than stresses in beam 3, case 2. Friction 

1/2 percent. less than 1/2 percent. Third approximation. 
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Figure 49.- Examples for idealization of cross sections, 
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Figure 44.- Comparisons between calculated and experimental stresses 
in beam of reference 5 for tip load. 
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Figure 45.- Comparisons between calculated and experimental stresses 
in beam of reference 5 for distributed load. 
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Figure 47.- Comparisons between calculated and experimental stresses for 
beam 1, with small cut-outs. 
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Figure 48.- Comparisons between calculated and experimental stresses for 
beam 1, with large cut-outs. 
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Figure 50.- Beam used for numerical examples. 
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Figure 51." Be*m used toy Ibner-lblle* for numerical example 
(from reference 13). 9/1 - 0.385. Dimensions 
are in centimeter unite. 
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Figure 53.- Ebner-i.611er beam, load caee 1, analyzed by different aetbode. 
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